CLASSIFICATION OF GLOBAL AND BLOW-UP 
SIGN-CHANGING SOLUTIONS OF A SEMILINEAR HEAT 
EQUATION IN THE SUBCRITICAL FUJITA RANGE II. 
HIGHER-ORDER DIFFUSION 

V.A. GALAKTIONOV, E. MITIDIERI, AND S.I. POHOZAEV 

Abstract. A detailed study of two classes of oscillatory global (and blow-up) solutions 
was began in |20| for the semilinear heat equation in the subcritical Fujita range 

(0.1) ut = Au + \u\P-^u in x R+ for 1< p < po = 1 + 

with bounded integrable initial data u(x, 0) = uq(x). This study is continued and ex- 
tended here for the 2mth-order heat equation, for m > 2, with non- monotone nonlinearity 

(0.2) ut = -(-A)™u + \u\P in x M+, in the range 1 < p < po = 1 + 

with the same initial data uq. The fourth order bi- harmonic case to = 2 is studied in 
greater detail. The blow-up Fujita-type result for (10.21) now reads as follows: blow-up 
occurs for any initial data uq with positive first Fourier coefficient: 

/ uo{x) dx > 0, 

i.e., as for (jO.ip . any such arbitrarily small initial function uo{x) leads to blow-up. The 
construction of two countable families of global sign changing solutions is performed on 
the basis of bifurcation/branching analysis and a further analytic- numerical study. In 
particular, a countable sequence of bifurcation points of similarity solutions is obtained: 

pz-l + ffj, / -0,1,2,.... 



1. INTRODUCTION: HIGHER-ORDER SEMILINEAR HEAT EQUATIONS, BLOW-UP, FUJITA 
EXPONENT, AND GLOBAL OSCILLATORY SOLUTIONS 

This paper generalizes to higher-order (poly-harmonic) diffusion operators the study of 
[20] of the semihnear second-order heat equation in the subcritical Fujita range 

(1.1) Ut = Au + lulP-'^u in X M+, where 1 < p < po = I + j^, 
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with bounded integrable initial data u{x,0) = uq{x). Necessary key references, results, 
and our main motivation of the study of (11.11) and further related models can be found 
in [201 § !]• 

1.1. Higher-order semilinear heat equation: blow-up Pujita-like result. Thus, we 
intend to extend some of key results of [20j to semilinear parabolic equations with higher- 
order diffusion operators. Such models are steadily becoming more and more popular 
in various applications and in general PDE theory. Namely, we consider the 2mth-order 
heat equation, for m > 2, in the subcritical Fujita range: 

(1.2) = -(-A)™u+ in X M+ for 1 < p < po = 1 + ^, 
with sufficiently smooth, bounded, and integrable initial data, 

(1.3) u{x,0)=Uoix) in M^. 

The choice of the non-monotone nonlinearity \u\^ (a source term) in (II. 2p is associated 
with the necessity of having a standard sounding Fujita- type blow-up result. Namely, it is 
known [8] that blow-up occurs for (II. 2p for any solutions with initial data having positive 
first Fourier coefficient (see [9] for further details and |21] for an alternative proof): 

(1.4) /jgjv uo{x) dx > 0, 

i.e., any, even arbitrarily small, such data lead to blow-up. A similar result for positive 
solutions of (II. ip was well-known since Fujita work in 1966; see [201 § 1] fo^' ^ survey. 

1.2. Results and layout of the paper. In Sections [SHU we perform construction of 
countable sets of global sign changing solutions on the basis of bifurcation/branching 
analysis, as well as of a centre-stable manifold one. Here, we apply spectral theory of 
related non-self-adjoint 2mth-order operators [9J, which is available for any m = 2,3, ... . 
As for (II. ip . i.e., for m = 1, this gives a similar sequence of critical bifurcation exponents: 

(1.5) p, = l + |^, / = 0,1,2,.... 

References and some results for analogous global similarity solutions of a different 
higher-order reaction-diffusion PDE with a standard monotone nonlinearity as in (II. ip . 

(1.6) ut = -{-A)"'u + \u\P-\ (m>2) 

can be found in [T?]. It is remarkable (and rather surprising for us) that the bifurcation- 
branching phenomena therein for (11.61) are entirely different from those for the present 
equation (11.20 . which turn out also to be more complicated, with various standard and 
non-standard bifurcation phenomena. 

It is worth mentioning here that our study also directly concerns blow-up solutions of 
(II. ip : we claim that, under the conditions that our two classes of its global oscillatory 
solutions are evolutionary complete (see |20i § 7] for a precise statement and some results 
for (II. ip ). all other solutions of (11.20 must blow-up in finite time. Then this describes a 
much wider class of blow-up solutions, and actually says that almost all (with a. a. defined 
in a natural way) solutions of (ll.2p in the subcritical Fujita range blow-up in finite time. 
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2. Global similarity solutions and ^-bifurcation branches 

In what follows, we use a general scheme and "ideology" of the study in |2(lj of the 
second-order semilinear equation fll.ip . Therefore, omitting some obvious details, we now 
more briefly start to describe which results on global solutions can be extended to the 
2mth-order react ion- diffusion equation (11 ■2p in the subcritical range. 

2.1. First family of global patterns: similarity solutions. As usual, for the higher- 
order model (11.21) with m > 2, we first study the existence and multiplicity of the standard 
global (i.e., well defined for all t > 0) similarity solutions of the form 

(2.1) us{x,t) =r^f{y), y = x/t^. 

This leads to the semilinear elliptic problem for the rescaled similarity profile /: 




BJ+\fr^-{-Arf + ^yVf + ^,f + \f\^ = in M 
f{y) decays exponentially fast as \y\ ^ oo. 



For m = 1, this problem admits a variational setting in a weighted metric of L^(]R^), 

where p = e'^'^/'*. This positive fact was heavily used in [20], where category/fibering 
techniques allowed us to detect a countable number of solutions and bifurcation branches. 

However, for any m > 2, (12. 2p is not variational in any weighted space; cf. reasons 
for that and a similar negative result in [231 § 7]. So, those power tools of potential 
operator theory in principle cannot be applied for (12. 2p . with any m > 2. 

Moreover, unlike the previous study of (II. ip in |i20j, we cannot use standard variational 
results on bifurcation from eigenvalues of arbitrary multiplicity (for our purposes, the 
results for odd multiplicity [71 p. 381, 401] concerning local and global continuation of 
branches are sufficient). We also do not have global multiplicity results via Lusternik- 
Schnirel'man (L-S, for short) and fibering theory. As usual, higher-order semilinear elliptic 
equations such as (12. 2p . or even the corresponding ODEs for radially symmetric profiles /, 
become principally different and more difficult than their second-order variational coun- 
terparts. We again refer to Pl § 6.7C] for general results on bifurcation diagrams, and to 
[Tl [29j for more detailed results for related 2mth-order ODEs in ID. These results do not 
apply directly but can be used for a better understanding of global bifurcation diagrams 
of similarity patterns f{y)- 

Thus, as in [17] for a quite similar looking equation (II. 6p . for global continuation of 
branches, we have to rely more heavily on numerical methods, and this is an unavoid- 
able feature of such a study of nonlinear higher-order equations. Surprisingly, we detect 
completely different local and global properties of p-branches in contrast with those in 
[T7] for the equation (II. 6p . which therefore are not so definitely attached to variational, 
monotone, or order-preserving (i.e., via the Maximum Principle) features of these difficult 
global similarity problems studied since the 1980s. 



2.2. Fundamental solution and Hermitian spectral theory. We begin with the nec- 
essary fundamental solution 6(x, t) of the corresponding linear parabolic (poly-harmonic) 
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equation 

(2.3) ut = -{-A)"'u in x M+, 
which takes the standard similarity form 

(2.4) h{x,t) = t^^^F{y), where y = x/t'^. 

The rescaled kernel F is then the unique radial solution of the elliptic equation 

(2.5) BF = -{-A)"'F + 2^ y ■ VF + ^ F = in M^, with J F = 1. 

The rescaled kernel F{\y\) is oscillatory as |y| — t- oo and satisfies the estimate [TUl E] 

(2.6) \F{y)\ < D e-"^!^!" in M^, where a = e (1,2), 

for some positive constants D and d depending on m and N. The linear operator Bi in 
equation (12. 2p is connected with the operator (12. 5 p for the rescaled kernel F in (12. 4p by 

(2.7) Bi = B + ciJ, where ci = ^-^^^ and Po = 1 + 

In view of (12. 7p . in order to study the similarity solutions, we need the spectral proper- 
ties of B and of the corresponding adjoint operator B*. Both are considered in weighted 
L^-spaces with the weight functions induced by the exponential estimate of the rescaled 
kernel (12. 6p . For m > 2, we consider B in the weighted space L^(R^) with the exponen- 
tially growing weight function 

(2.8) p{y) = e'^l^l" > in M^, 

where a G (0, 2d) is any fixed constant and d is as in (12. 6p . We ascribe to B the domain 
i7^"*(M^) being a Hilbert space with the norm 

2m 



\V\ 



|2 



Jpiy)Z\D'viy)\'dy, 



k=0 

r2m 



induced by the corresponding inner product. Then i/;™ C C L\ The spectral 
properties B follows [9]: 

Lemma 2.1. (i) B : H^"^ L'^^ is a bounded linear operator with the real point spectrum 

(2.9) a{B) = {\i = -£, / = 0,1,2,...}. 

The eigenvalues Xi have finite multiplicities with eigenf unctions 

(2.10) My) = ^ D^Fiy), for any = /. 

(ii) The set $ = {'?/'/3}|/3|>o is complete and the resolvent (B — XI)^^ is compact in L^. 

By Lemma [2.11 the centre and stable subspaces of B are given by E'^ = SpanjT/'o = -^} 
and E' = Span{?/';3, > 0}. 

Consider next the adjoint (in the dual L^-metric) operator 

(2.11) B* = -(-A)'"-^yV. 
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For m > 2, we treat B* in L^* with the exponentially decaying weight function 

(2.12) p*(^) = _L_^e-l^l">0. 

Lemma 2.2. (i) B* : H'p^ — )■ L^* is a hounded linear operator with the same spectrum 
02. 9p as B. Eigenf unctions ip*p{y) with = / are Ith-order generalized Hermite polyno- 
mials given by 

m/2m] 

(2.13) r,{y) = :}p^, E l!(-A)-V 



(ii) The set $* = {'?/'«}|/3|>o ^-5 complete and resolvent (B* — XI) ^ compact in L'] 



It follows that the orthonormality condition holds 

(2.14) (V'/3,^;) = 5/3^, 

where (■, ■) denotes the standard L^(R^) inner product and 5/?^ is Kronecker's delta. 

Using fl2.14p . we introduce the subspaces of eigenfunction expansions and begin with 
the operator B. We denote by the subspace of eigenfunction expansions f = E C/?^'/? 
with coefficients c^ = {v,ip*) defined as the closure of the finite sums {E|^|<m '^^^/^i 
the norm of L^. Similarly, for the adjoint operator B*, we define the subspace L^* C L^, . 
Note that since the operators are not self-adjoint and the eigenfunction subsets are not 
orthonormal, in general, these subspaces can be different from and L^* , and the equality 
is guaranteed in the self-adjoint case m = 1, a = j only. 

2.3. Existence of similarity profiles close to transcritical bifurcations. Consider 
the elliptic problem (12. 2p . Using the above Hermitian spectral analysis of the operator 
pair {B,B*}, we formulate the bifurcation problems, which guarantee the existence of a 
similarity solution in a neighbourhood of bifurcation points. In fact, our consideration is 
quite similar to that for the second-order case in [20j, so we may omit some details. Since 
p < Pq, our analysis is performed in the subcritical Sobolev range: 

(2.15) l<p<ps = =^ iJ™(M^) c LP+i(M^) compactly. 

Taking p close to the critical values, as defined in (II. 5p . we look for small solutions of 
(12. 2p . At p = pi, the linear operator Bi has a nontrivial kernel, hence: 

Proposition 2.1. Let for an I > 0, the eigenvalue = operator (12. 5p is of odd 

multiplicity. Then the critical exponent (II. 5p is a bifurcation point for the problem ( 12. 2p . 

Proof. Consider in our equation written as 

(2.16) B/ = -(1 + ci)/ - l/r, where B = Bi - (1 + ci)J = B - J. 

It follows that the spectrum o"(B) = {—1 — ^} consists of strictly negative eigenvalues. 

The inverse operator B~^ is known to be compact, [9l Prop. 2.4]. Therefore, in the 
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corresponding integral equation 

(2.17) / = A(/) ^ -(1 + ci)B-V - B-\\fr), 

the right-hand side is a compact Hammerstein operator; see [26l Ch. V] and apphcations 
in [3l[T7l[23]. In view of the known spectral properties of B~^, bifurcations in the problem 
(12.170 occur if the derivative A'(0) = — (l+Ci)B^^ has the eigenvalue 1 of odd multiplicity, 
[271 126]. Since ct(A'(0)) = {(1 + ci)/(l + ^)}, we obtain the critical values (O]). By 
construction, the solutions of f l2.17p for p ^ pi are small in and, as can be seen from the 
properties of the inverse operator, / is small in the domain Hp"^ of B. Since the weight 
(12. 8p is a monotone growing function as \y\ — )■ oo, using the known asymptotic properties 
of solutions of (12. 2p . / G Hp"^ is a uniformly bounded, continuous function (for N < 2m, 
this directly follows from Sobolev's embedding theorem). □ 

Thus, Z = is always a bifurcation point since Aq = is simple. In general, for 
/ = 1,2, ... the odd multiplicity occurs depending on the dimension A^. For instance, for 
1 = 1, the multiplicity is A^, and, for / = 2, it is i^I^^±ll. in the case of even multiplicity 
of A;, an extra analysis is necessary to guarantee that a bifurcation occurs |27j using the 
rotation 7/ of the vector field corresponding to the nonlinear term in (I2.17P on the unit 
sphere in the eigenspace = Spanj?/'/?, = 1} (if 7^ 7^ 1, then bifurcation occurs). We 
do not perform this study here and note that the non-degeneracy of this vector field is 
not straightforward; see related comments below. It is crucial that, for main applications, 
for = 1 and for the radial setting in M^, the eigenvalues (12. 9 p are simple and (II. 5p 
are always bifurcation points. Unlike Proposition 3.2 for m = 1 [20, § 3], we have the 
following result describing the local behaviour of bifurcation branches occurring in the 
main applications, see [26] and [271 Ch. 8]. Unlike the case m = 1, some bifurcations 
become trans critical. 



Proposition 2.2. Let Xi be a simple eigenvalue o/B with eigenfunction ipi, and let 

(2.18) «:, = (|^,|P,^;)^0. 

Then the p-bifurcation branch crosses transversely the p-axis at p = pi. 

We next describe the behaviour of solutions for p ^ pi and apply the classical Lyapunov- 
Schmidt method, [271 Ch. 8], to equation (I2.17P with the operator A that is differentiable 
at 0. Since, under the assumptions of Proposition 12.21 the kernel Eo = kerA'(O) = 
Span{'i/'/} is one- dimensional, denoting by Ei the complementary (orthogonal to ipi) in- 
variant subspace, we set 

(2.19) f = Fo + Fi, where Fq = eitPi e Eo and Fi = Zk^^Sk^k e E^. 

Let Pq and Pi, Pq + Pi = I, be projections onto Eq and Ei respectively. Projecting (I2.17P 
onto Eq yields 

/o on\ , /-D-l/l^lPA „;.*\ 1 l+ci (N+l)s „ „ „ 



By bifurcation theory (see [271 P- 355] or [71 p. 383], where A'(0) is Fredholm of index 
zero), Fi = o{ei) as ei — )■ 0, so that ei is calculated from ( 12.20p as: 

(2.21) ^^el = -\El\P{B-^\^Pl\P,^P|)+oief) =^ le.r^^, = q(p - Pz)[l + o(l)], 
where q = '■^^^ • We have used the following calculation: 

(B-i|^,p,^;) = (i^d^(B*)-vr) = -li- 

Recall the identity (B-^)* = {B*)-\ 

It follows from the algebraic equation in f l2.2ip that the bifurcations are transcritical 
provided that ki ^ 0, while the sign of ki determines how the branches cross the p-axis. 
For / = 0, it is easy to see that > 0, since ipo = F and ipQ = 1, so that 

(2.22) Ko = {\H'',ro)=I\F\^>0 {p = Po). 

Moreover, for pq ^ 1+, by the definition of F in f l2.5p . we have that ^ 1. The 
positivity or negativity of the scalar product f l2.18p for / > 1 and arbitrary p > 1 is not 
straightforward, and we should rely on a delicate numerical evidence; see [T7]. It turns 
out that Ki can be both positive and negative for different / > 1. 

Let us note the following principal difference in comparison with the case of monotone 
nonlinearity +\f\^~^f studied before for m = 1, pOj. It turns out that 

(2.23) Ki = for odd / = 1,3,5,..., 

since in f l2.18p . \il'i{y)\^ is even, while the polynomial ipiiu) is odd. This means the 
following: 

(2.24) P = Pi for odd / > 1 are not "standard" bifurcation points. 

The corresponding "non-standard" bifurcation phenomenon will be discussed shortly. On 
the other hand, for the standard monotone nonlinearity as in ( II. 6p . all critical exponents 
{pi} are pitchfork bifurcations, [T7] . 

Thus, under the assumption ki ^ on the coefficients fl2.18p . we obtain a countable 
sequence of bifurcation points (11. 5p satisfying j9; — )■ 1+ as / — )■ oo, with typical transcritical 
bifurcation branches appearing in a neighbourhood. The behaviour of solutions in H^"^ 
and uniformly in M^, for p < ps, takes the form 

(2.25) fi{y) = \ci{p-pi)\~ sign {ci{p-pi)){'ipi{y) + o{l)) as p pi. 

Instability of all these local branches of similarity profiles is studied similar to the case 
m = 1 in ^ § 2]; see also [IT]. 

2.4. Lyapunov— Schmidt branching equation in the general multiple case: non- 
radial patterns. Let now = have multiplicity M = M{1) > 1 given by the 
binomial coefficient 

(2.26) M{1) = dim W%B - Xil) = = fj^, so that 



(2.27) 



Eo = ker(B - A//) = Span{^n, ...,^/m}- 
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Then, looking for a solution 

(2.28) f = fo + /i, with /o = eiipn + ... + Em^im, where fi-LEo, 

and substituting into the equation fl2.17p . multiplying by ipj^, and denoting, as usual, 
s = p — Pi, < |s| <^ 1, we obtain the following generating system of M algebraic 
equations: 

(2.29) e, = j^^ J\ei^n + ...+eMiJiM\''i^u = D,{eu...,eM), ^ = 1,2,...,M. 

Here p = pi. Denoting x = (ei, ...,£:m)"^ € K*^, the system fl2.29p is written as a fixed 
point problem for the given nonlinear operator D = {Di{x), Dm{x))^ , 

(2.30) X = B{x) in . 

In the second-order case m = 1 |20, § 3], the system (12.301) was variational, that allowed 
us to get a multiplicity result. In view of the dual metric in (I2.14p . for any m > 2, 
the algebraic system (12.291) is not variational, so the multiplicity of admissible solutions 
remains an open problem. 

Global extension of the above local p-bifurcation branches is performed by classic theory 
of nonlinear compact operators, [TK^SIEZ]- However, since the problem is not variational, 
nothing prevents existence of closed p sub-branches or appearances of turning, saddle- 
node bifurcations (we will show that this can actually happens), so that the total number 
and structure of solutions for any p G (l,Po] remain a difficult problem. We will then 
inevitably should rely on careful numerics. 

2.5. "Non-standard" pitchfork bifurcations for ki = 0. Without loss of generality, 
we consider the simplest case I = 1, N = 1, m = 2 (then pi = 1 + m = 3 by (II. 5p ). where, 
from (I2.18P and (I2.10p . (I2.13p . it is clear that ki vanishes: 

(2.31) /ti = {\H''\ri) = /m \F'{y)\hdy = 0. 

Next, unlike the standard approximation (I2.19P close to p = 3, we now use an improved 
one given by the expansion on the 2D invariant subspace E12 = Span{-?/;i, ■?/'2} (this choice 
will be explained below): 

(2.32) / = Fi2 + F3, where F12 = ei^i + £2^2 and ^3X^12, 

with the scalar parameters ei, £2 to be determined. For simplicity, we next use the 
differential version of the integral equation (I2.17P for / = 1: 

(2.33) (B-Ai/)/-f /= + where s = p - p, 

and where we omit the 0(s^)-term. Substituting (I2.32p into (I2.33P and projecting onto 
corresponding one-dimensional subspaces, quite similar to the system (12. 29 p . we obtain 
the following asymptotic system of two algebraic equations: 



(2.34) 



'1^2 = 1^2 - / klV"! +^2V'2|^^2 + 
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where we put -^2 = — 2m ~ ~I (^^^ I = m = 2) and where we have omitted higher- 
order terms associated with the orthogonal in (12.321) and via replacing p hj pi = 3 
in the integrals on the right-hand sides. Then, the second equation, as s — > 0, gives the 
dependence of 62 on the leading expansion coefficient ei on £'12: 

(2.35) £2 = 2|£:ip/ii2 + where /ii2 = / lV'iPV'2 7^ and ^'2 = 752/^- 

It is crucial that, unlike in (]2.3ip . the coefficient fiu is given by the integral of some even 
function, so that, now, the assumption /ii2 7^ is not that restrictive and can be quite 
reliably checked numerically. 

Next, the ffist equation in (12.341) . after a Taylor expansion in the integral, by using that 
€2 = o{ei) as s — 7- 0, provides us with the necessary bifurcation scalar equation on ei, 

(2.36) fei = J \eiipi fipl + 3el sign 61621^12 + ■-, with U12 = J ijf{signipi) i;2'ipl ^0, 

where, again, in 1^12, we face an even function in the integral. Since the ffist coefficient 
vanishes by (12.231) for I = 1, J lipifipl = 0, using the dependence (I2.35p . we obtain 

(2.37) \ei\^ = Ci2S + where 012 = 24^^. 

It follows that we thus deal with a pitchfork bifurcation at p = pi = 3, which is subcritical 
if C12 < and supercritical if C12 > 0. 

Overall, the bifurcation branches take the following form: for, e.g., C12 > 0, 

(2.38) f{y) = ±[c,2{p-S)]-sMy) + 2fiu[ci2{p-S)]lMy) + - as p^3+. 

We will reveal this kind of bifurcation numerically in Section HI Note that this non- 
standard bifurcation branch near p = 3 is more "steep", ~ 0((p — 3)5)), than the standard 
one in (I2.25p . which, for / = 1, is of the order ~ 0{y/p — 3). 

One can see that a similar bifurcation scenario, under the vanishing assumption (I2.23p . 
can be developed by using other invariant subspaces rather than that in (12. 32 p . The crucial 
conditions then remain the same: the corresponding coefficients /x., in (I2.35P and u in 
(I2.36P must be non-zero, which is possible by mixing even and odd eigenfunctions in the 
subspace, depending on the multiindices chosen. This has an interesting and surprising 
consequence: 

(2.39) (I2.23P : there can be more than one bifurcation branch, even for ID eigenspace. 
In Section m we will observe this numerically for the one-dimensional eigenspace. 

2.6. Transversality of intersections of subspaces. This was a permanent subject of 
an intense study for nonlinear second-order parabolic equations; see related key references 
and further comments in [201 § 6.2]. We briefly recall these important results below. 
Namely, this problem was completely solved rather recently for a scalar reaction-diffusion 
equation on a circle of the form 

(2.40) ut = A{u) = + g{x, u, u^), x G 5*^ = M/27rZ, 

where the nonlinearity g{-) satisfies necessary conditions for existence of global classical 
bounded solutions for arbitrary bounded smooth initial data. Then, if / is a hyperbolic 
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equilibrium of A, A(/) = 0, known to be generic (or a rotating wave), then the global 
stable and unstable subspaces of A'(/) span the whole functional space = H'^°'{S^), 
a E (|) 1)5 where the global semiflow is naturally defined, i.e., 

(2.41) W'iA'if)) © W{A'{f)) = X", 

so that these subspaces intersect transversely. It is crucial that such a complete analysis 
can be performed in ID only, since it is based on Sturmian zero set arguments (see [16] 
for main references and various extensions of these fundamental ideas), so, in principle, 
cannot be extended to equation in M^. We refer to most recent papers O [131 [21], where 
earlier key references and most advances results on the transversality and connecting 
orbits can be found. 

We perform our transversality analysis for p close to the bifurcation points p ^ pi in 
fll.Sp by using bifurcation theory from Section 12.31 

Proposition 2.3. Fix, for a given p ^ Pi, p ^ Pi, a hyperbolic equilibrium fj^, with a 
= /, of the operator A in (12. 2p . 

(2.42) A(/) = -(-A)-/ + A. ^ . V/ + ^ / + 
Then the transversality conclusion holds: 

(2.43) W'{A:{fp)) © W-{A!{fp)) = //^-(mA^). 

Proof. It follows from (I2.42p and the expansion fl2.25p that, for p = pi+e, with < |£| ^ 1, 
A'(/^) = -(-A)- + J- 2/ . V + ^ J + sign fp I 

= (B - \il)+p\ci\ \e\\^i\^-hign{cie^p) + ... (/ = |/3|) . 

Therefore, for p = pi, the following analogy of f l2.43p is valid: 

(2.45) A'ifp) = B-Xil =^ W'{B-XiI)®W''{B-XiI)®W''{B-XiI) = i/|'"(M^), 

and diml^'^(B — A//) is equal to the algebraic multiplicity f l2.26p of A/ = By the as- 

sumption of the hyperbolicity of and in view of small perturbations (see, e.g., [ll[2n]) of 
all the eigenfunctions of A'(/^) for any \6\ <^ 1, e 0, which remain complete and closed 
as for p = pi, we arrive at f l2.43p . Recall that, since by f l2.44p . A\fp), with eigenfunction 
l'?/'^}, is a small perturbation of B — A^/ (with eigenfunctions {^Jp}) and, in addition, 
the perturbation is exponentially small as ?/ — )■ 00, the "perturbed" eigenfunctions ippiy) 
remain a small perturbation of the known ipp{y) in any bounded ball, and sharply approx- 
imate those as ?/ — > 00. Therefore, close to p = pi, there is no doubt that the well-known 
condition of completeness/closure of {ipis} (the so-called property of stability of the basis) 
is, indeed, valid: 

Thus, close to any bifurcation point p = pi, we precisely know both the dimensions of 
the unstable subspace of AH^fp) of any hyperbolic equilibrium fp (and, sometimes, we can 
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prove the latter) and the corresponding eigenfunctions {ipp}- 

(2.46) by continuity, for all p ~ : -^^ ~ ^-^i = im i^p ^ i^p^ 

where convergence of eigenfunctions as p — )■ holds in iY^™ and uniformly in M^. 

Furthermore, moving along the given bifurcation p-branch, the transversality persists 
until a saddle-node bifurcation appears, when a centre subspace for A'(/^) occurs, and 
hence fl2.43p does not apply. If such a "turning" point of a given p-branch does not appear 
(but sometimes it does; see Section H] below) , the transversality persists globally in p. 



3. Numerical results: extension of even p-branches 

Thus, the above bifurcation analysis establishes existence of a countable set of transcrit- 
ical p-bifur cat ions aX p = pi for even /.As we have mentioned, since (12. 2p is not variational 
for m > 2, we do not have any chance to use power tools of category-genus- fibering the- 
ory in order to guarantee nonlocal extensions of p-branches of similarity profiles f{y)- 
However, as is well-known from compact nonlinear integral operator theory |71 |26l 127] , 
these branches are always extensible, but can end up at other bifurcation points, so their 
global extension for all p > pi is not straightforward. Actually, we show that precisely 
this happens for m = 2 in ID. 

3.1. Preliminaries for m = 2: well-posed shooting of even profiles. We first 
concentrate on the simplest fourth-order case: 

(3.1) N = 1 and m = 2, so that po = 1 + ^ = 5, 

in order to exhibit typical difficult and surprising behaviours of global p-branches of the 
first similarity profile fo{y), which bifurcates from the first critical exponent po = 5 in 
(13. ip . We also compare /o in dimensions A = 1, 2, 3, and 4. For convenience, we will 
denote by fi{y) the profiles that bifurcate at the corresponding critical pi and hence, by 
(I2.25p . "inherit" the nodal set structure of the eigenfunction ipi{y) in (I2.10p for A = 1. 
In the case (13. ip . the problem (12. 2 p becomes an ODE one: 

^32) 1 A(/)^-/(^) + iy/' + ^/+|/|^' = for y>0, 

yf{y) decays exponentially fast as ?/ — )■ ±oo. 
We first easily prove the following result, somehow confirming our bifurcation analysis: 

Proposition 3.1. (i) In the critical case p = 5, the only solution of ( 13. 2p is f = 0; and 
(ii) The total mass of solutions of (13. 2p satisfies 

(3.3) J f < for p < Po = 5 and J f > for p > 5. 
Proof. Integrating the ODE (13. 2p over M yields the following identity: 

(3.4) I\f\'=^J'f- □ 
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Remark on bifurcation analysis. Firstly, according to the bi-orthogonality (I2.14p . 

(3.5) Jipi = for all / = 1,2,3,..., and J = J F = 1, 

so we see that (12.251) somehow "contradicts" (13. 3p . However, there is no any controversy 
here: indeed, (12.19^ assumes, in ID, the following expansion: 

(3.6) f = eiipi + eoipo + ... , 

where we keep the only eigenfunction tpQ with the unit non-zero mass. Then, the identity 
(13. 4p is perfectly valid provided that 

(3.7) eo = '-^\ei\Pj\mi + o{l))=o{e,), 

so that this small correction in (13. 6p allows one to keep the necessary non-zero mass on 
any even p-bifurcation branch. 

As usual, for the even profile /o (and for /2, f^,...), since the ODE (13. 2 p is invariant under 
the symmetry refiection y — two symmetry conditions at the origin are imposed, 

(3.8) /'(O) = r(0) = (then f{-y) ^ f{y)). 

Let us first reveal a natural "geometric" origin of existence of various solutions of the 
problem (13. 2p . (13. 8p . This is important for the present non- variational problem, where 
we do not have other standard techniques of its global analysis. It is easy to see that the 
ODE in (13. 2p admits 2D bundle of proper exponential asymptotics as y ^ +oo: 

(3.9) f{y) ~ e-'^y''' [d cos yl) + sin , ao = 3 ■ 2-f , 

where Ci^2 G are arbitrary constants. Obviously, (13. 2p also admits a lot of solutions 
with much slower algebraic decay, 

(3.10) /(2/)~Coi/"^ as y^+oo, CqER, Co ^ 0, 

but these should be excluded from the consideration, so we always must take Cq = 0. 

These two parameters Ci^2 in (13. 9 p are used to satisfy (to "shoot") also two condi- 
tions at the origin (13. Sp . Overall, this looks like a well-posed ("2—2", i.e., not over- and 
under-determined) geometric shooting problem, but indeed extra difficult "oscillatory" 
properties of the ODE involved are necessary to guarantee a proper mathematical con- 
clusion on existence of solutions and their multiplicity (in fact, an infinite number of 
those). This will be done with the help of numerical methods, and, as was mentioned, the 
final conclusions are striking different from those obtained in [T71 [181 HH] for monotone 
nonlinearities. 

Thus, we arrive at a well posed "2 — 2" shooting problem: denoting by / = f{y; Ci, C2) 
solutions having the asymptotic behaviour (13. 9p (note that such solutions can blow-up at 
finite ?/o > 0, but we are interested in those with yoiCi, C2) < 0; see below), by (13. 8p . an 
algebraic system of two equations with two unknowns occurs: 

.o^^^ r/'(0;Ci,C2) = o, 

> \f"(0;Ci,C2)=0. 
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Proposition 3.2. For any even integers p = 2,4,6, the system (13. lip admits not more 
than a countable set of solutions. 

Proof. For such p's, the ODE fl3.2p has an analytic nonhnearity, so by classic ODE theory 
[6l Ch. I], both functions in (13. lip are also analytic, whence the result. □ 

We expect that a similar result is true for arbitrary p > 1, but a proof of an analytic 
dependence on parameter^], is expected to be very technical. 

3.2. The first symmetric profile fo{y). For solving our problem (13. 2p . we use the 
bvp4c solver of the MatLab with the enhanced accuracy and tolerances in the range 

(3.12) 10"^ - 10"^^ 

and always, with a proper choice of initial approximations (data), observed fast conver- 
gence and did not need more than 2000-8000 points, so that each computation usually 
took from 15 seconds to a few minutes. 

We begin with Figure [1] presenting a general view of the similarity profile fo{y) for 
various p above the critical exponent po = 5. It is clearly seen that fo{y) is oscillatory for 
large y, but definitely has a dominated "positive hump" on y E (0,3.4), so that overall 
(cf. O) 

//o>0. 

However, by (11. 4p . this does not imply blow-up of the corresponding similarity solution 
u{x,t), since this happens in the supercritical range p > po = 5, when, in particular, all 
sufficiently small solutions are known to be global in time. 

Figure [2] shows the dependence of the radial pattern /o = fo{\y\) on dimensions = 1, 
2, 3, 4. All the profiles look similar and their L°°-norm, /o(0), increases with A^. However, 
the location of the "positive hump" of each fo^y) remains practically unchanged, as well 
as the location of the first "nonlinear transversal zero" , ?/o ~ 4 always; see more below. 

3.3. p-branches and further even profiles. More delicate results are shown in Figure 
El where we present the global po-branch, initiated at p = 5+ and extended up to p = 200. 
In particular, this shows that 

(3.13) ll/olloo = /o(0) ^ 1+ as p^+oo, 

an asymptotic phenomenon with a possible difficult logarithmically perturbed behaviour 
that was discovered and studied in [T71 § 5] for another model (13. 2 p with the monotone 
nonhnearity |/|^~^/. The deformation of the profile /o on the same interval p G [5.01, 200] 
is shown in Figure HI again confirming (I3.13P . 

We next study the extension of the po-branch for p < po = 5. The transition through 
the first transcritical bifurcation at p = po = 5 is explained in Figure |5l which shows 
a clear spatial similarity of foiv) ~ ^^oiu) = ^F{y) along both limits p — )■ Po =5^, 
according to (I2.25P for / = 0. 

^As is well known, dependence on parameters in such ODE problems can be much better than the 
smoothness of coefficients involved. A classic example is: for elliptic operators with just measurable 
coefficients, the resolvent is often a meromorphic function of the spectral parameter A € C. 
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global similarity profile ^^(y): N=1 , m=2, p s [5.01 , 6.4] (Pq=5) 



0.7 



f(y) 




2 4 6 8 10 

y 

Figure 1. The first profile fo{y) of ([SJ]) for iV = 1, m = 2 and p G [5.01, 6.4]. 




y 

Figure 2. The first radially symmetric solution fo{y) of (j2.2p for m = 2, p = 6, 
and dimensions = 1, 2, 3, 4. 
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N=1 , m=2; transition of f.(y) ttirough p =5 




y 

Figure 5. Transition of fo{y) to /2(y) of ([32]) for iV = 1, m = 2 for p 5=^. 

The global p2-brancli, which is an extension of the positive po-one in Figure HI is shown 
in Figure [6] while the corresponding deformation of /'s in Figure [71 It turns out that it 
ends up at the next (even) bifurcation point 

(3.14) p = p2 = 1+1^ = 1 = 2.3333..., 

so that the branch is expected to be continued for p < p2 = in a. "positive" way, etc. 

To justify such transcritical bifurcations at p = pi for even Z > 2, in Figure IH we present 
a transition through p = P2 = | • Similarly, in Figure [9l we show transition from fi^y) for 
p = 1.85 > p4 = 1 + I = 1-8 to faiy) for p = 1.75 < p^. 

Thus, according to the results given above, we expect that there is a continuous defor- 
mation along each connected branch of /o into /2, /2 into f^, into /g, etc., i.e., there 
exists a unique global p-branch of even similarity profiles. Hence, we observe that all 
connected branches have similar shapes with always two bifurcation points involved: the 
right-hand end point p = p^k and the left-hand end one ^4^+2- 

3.4. "Approximate" Sturmian zero property. Let us comment on the "Sturmian 
property" of similarity profiles {fiiy)}- Figures [5] and [8] indicate that, regardless the 
oscillatory exponential tails, each profile fi{y) has a clear "approximate" ("nonlinear") 
Sturmian structure and exhibits I + 1 dominant extrema (meaning / "transversal" zeros 
in between). In a rigorous mathematical sense, such properties are known to hold for 
the second-order problems. In [20], where (11. II) was studied, Sturmian properties were 
connected with the category of the functional subset for each fi, being the corresponding 
min-max critical point of the functional, since the category assumes using the reflection of 
the functions (■) 1— )■ — (■) and hence the nodal sets of fi{y) gets more and more complicated 
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N=1 , m=2: negative Pj^-branch of f^^ ^(y) for p e (P2,Pg)=(7/3,5) 

f(0) ° 

-0.1 
-0.2 



-0.3 
-0.4 
-0.5 



2 2.5 3 3.5 4 4.5 5 

P 

Figure 6. p2-branch of /2 for = 1, m = 2, extended for p G (p2)Po)- 

N=1 , m=2: deformation of f^ on the negative p~-branch 

0.2 I < > < < < < > > 1 




y 

Figure 7. p-deformation of /2 from Figure [H N = 1, m = 2, p £ {p2,Po)- 




N=1 , m=2: transition tinrougti P2=7/3=2.333... 
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P=2.4>P2 
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Figure 8. Transition of for V = 2.4 > p2 into /4 for p = 2.25 < where 
spatial shape of both mimics the eigenfunction ip2{y) according to (j2.25p . 



N=1 , m=2: transition tlirough p^=9/5=1 .8 
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Figure 9. Transition of /4 for p = 1.85 > ^34 into /g for p = 1.75 < p4, where 
spatial shape of both mimics the eigenfunction ■04 (y) according to (j2.25p . 

as / increase^. However, in the present non- variational case, we cannot use even those 
rather obscure issues, though numerical evidence clearly suggests that the approximate 
Sturmian nodal properties persist in both variational and non-potential problems. Hon- 
estly, we do not know, which "mathematical/functional structures" can be responsible 



There is no still a rigorous treatment of such zero-set phenomena. 
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for such a hugely stable Sturmian-like phenomena, and this remains a challenging open 
problem of nonlinear operator theory. 

For the higher-order equations, some extra mathematical reasons for Sturmian prop- 
erties to persist in an approximate fashion are discussed in [T7i § 4.4]. These can be 
attributed to the fact that the principal part of (13. 2 p contains the iteration of two positive 
operators 

-Dt = -{-Dl){-Dl), 

and for such pure higher-order operators Sturm's zero property is true [11]. Then the 
linear perturbations affect non-essential zeros in the exponential tails only. No rigor 
justification of such a conclusion is available still. 



4. Towards odd non-symmetric profiles and their ^-branches 

4.1. An auxiliary discussion. We begin with noting that, for such fourth-order non- 
linear ODEs f l3.2p with a clearly principal non-coercive operator, for any p > 1, we expect 
to have, at least, a countable set of different solutions (as in [T71 [TSl ITU]). 

However, in all our previous studies of higher-order elliptic ODE problems, [T71 [181 HH] , 
exactly half of such solutions were odd functions of y (in ID; in the radial symmetry, 
obviously, odd profiles are not admitted). In the present case, the profiles fi{y), fs{y),... , 
are not odd (anti-symmetric), since the ODE (13.20 . unlike that for (11.60 . does not admit 
the corresponding symmetry 

(4.1) f^-f, -y- 

Bearing in mind the multiplicity bifurcation results in Section 1531 o^ig concludes that other 
profiles must be non-symmetric in any odd sense. Notice that precisely that explains the 
bifurcation expansion (I2.38p . where the leading term is odd via ipi{y), while there exists 
always an even correction via ip2{y)- 

We first check some analytical issues concerning such unusual similarity profiles. Thus, 
we shoot from y = +oo using the same bundle as (13.90 . with the coefficients C12, 

(4.2) f{y) e-^y'^'[C t cos {^y I) sin {^yl)], where = 3-2-1. 

Evidently, most of such solutions / = f{y;Ci,C2) will blow-up at some finite yo = 
yo{C^, C2) according to the following asymptotics: as y ^ y^, 

.4.. /(^) = |/ini + o(l)) =^ /(y) = Co(y-yo)"^(l + o(l)), 

where C^o^"^ = ^(^ + l) ( ^ + 2) + 3). 
Therefore, in order to have a global profile, we have to require that 

(4.4) yoiCt,Ct) = -oo. 

Once we have got such a global solution defined for all ?/ G M, we then need to require 
that, ai y = —00, the algebraic decay component (I3.10p therein vanishes, i.e., 

(4.5) Co-(C+,C+) = 0. 
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Figure 10. A non-symmetric profile fsiy) for p = 2.7, well satisfying the 
identity (15:11) . 

We, thus, again arrive at a system of two algebraic equation ( 14. 4p . ( 14. 5p . where the 
result of Proposition 13.21 applies directly to guarantee that the total number of possible 
solutions is not more than countable. 

In Figure [T0| as a first typical example, we show a non-symmetric "dipole-like" profile, 
denoted by fsi^y) (see below why such a subscript) for p = 2.7, which with a sufficient 
accuracy ~ 10~^ satisfies the identity ( 13. 4p . It turned out that this identity can be used 
as a "blueprint" for checking the quality in some worse-converging cases. 

4.2. p-branch of f^: from pi = 3 to a saddle-node bifurcation. Starting from the 
profile /s for p = 2.7 in Figure [IDl we perform a continuation in the parameter p to get to 
the bifurcation origin of this p-branch. Not that surprisingly (cf. Section 12751) . we observe 
in Figure [TlTa) that the corresponding bifurcation branch, with certain accuracy, goes to 
the odd bifurcation point (11.51) . with 1 = 1: 

(4.6) = 1 + I = 3. 

In this calculation, we take the continuation step Ap = 10~^ (then the calculation takes 
a couple of hours), so, as seen, we cannot approach closely to this bifurcation point. To see 
approaching p = pi more clearly, we, in addition, took the continuation step Ap = 10~^ 
(the calculation of the full branch then took about 16 hours) and observe approaching 
Pi = 3 up to p = 2.9991, so that existence of such a "forbidden" earlier bifurcation is 
without any doubt. Figure fllT b) clearly shows that, as p — )■ = 3, the profile fsi^y) 
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(the bold dashed hne at p = 2.998)) takes a typical "dipole" behaviour governed by the 
second eigenf unction from (12.101) : 

(4.7) My) = -F'iy), 

where i^(||/|) is the even "bi-harmonic Gaussian" satisfying (12. 5p for m = 2. Therefore, 
this /sd/) ~ C{p)il}i{y), with some (unknown still) constant C{p), looks like a standard 
dipole for the Gaussian for m = 1 



(4.8) i;,{y) = ^ye 



-2/2/4 



but it has oscillatory tails and further dominant positive and negative humps. 

Extending this p-branch of /s for p < 2.7, we observe existence of a saddle-node bifur- 
cation at some p = Ps-n, where we obtain the estimate 



(4.9) 2.6148 < Ps-n < 2.6149 



by using again the step Ap = 10~^. The profiles /a close to Ps-n are shown in Figure | 

The lower bifurcation branch of /s is shown in Figure [131 The corresponding upper 
bifurcation branch is shown in Figure [HI while the corresponding deformation of /s is 
presented in Figure [T5l 

Quite surprisingly, the upper bifurcation branch in Figure [TH ends up at the previous 
bifurcation point po = 5! Comparing with Figure [H] for symmetric even profiles, we thus 
obtain two different bifurcation branches (of symmetric and non-symmetric) solutions 
originated at p = 5~. It is worth mentioning that the kernel of the linearized operator at 
p = 5 is one- dimensional. We still do not have a proper explanation of such a hard and 
unusual bifurcation phenomenon. However, the case / = is not degenerate by (I2.22p . 
so it seems one cannot create a bifurcation approach similar to that in Section 12.51 This 
remains an open problems 
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Figure 12. Non-symmetric profiles /^{y) close to the saddle-node bifurcation (|4.9p . 
N=1 , m=2: p-bifurcation diagram of yy), the lower branch; A p=1 0"'^ 
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Figure 13. The lower bifurcation branch of the non-symmetric profile fsiy) 
for p G [2.6149,2.7]. 
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N=1, m=2: p-bifurcation branch of yy); the upper branch for ps [Pg_^>5] 




Figure 14. The upper bifurcation branch of the non-symmetric profile fsiy) 
for p G [2.6149,5]. 

N=1 , m=2: deformation of fg(y) for p s [Pg_^,5], A p=0.005 




Figure 15. Deformation of the non-symmetric profile fsiy) for p G [2.6149,5]. 
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5. Second countable family: global linearized patterns 

5.1. Stable manifold patterns. This construction is similar to that for m = 1 and 
again relies on stable manifold and Hermitian spectral theory for the operator pair 
{B,B*} in Section fI72\ We perform the same scaling of a global solution u{x,t) of f ll.2p 
for t > 1, 



(5.1) u{x,t)=t p-'^v(y,T), y = x/t2,n^ r = lnt, where 

(5.2) Vr = A{v) = -{-^)"'v + ^yVv + ^^v+\v\P, so 

(5.3) A'(0) = B + ci/, ci = ^-^ = ^M>0 for p<po. 
Thus, A'(0) has the infinite-dimensional stable subspace: 

(5.4) E' = Span{ij,3 : X^ + ciKO, i.e., \(3\ > 2ci}. 

Using the above spectral properties of B [9], similar to [20l § 5], by invariant manifold 



theory for parabolic equations [28l Ch. 9], we arrive at the following (see also [T5l [17] ): 

Proposition 5.1. For any multiindex f3 satisfying \f3\ = I > 2ci, equation (15.21) admits 
global solutions with the behaviour, as t ^ +oo, 

(5.5) t;;3(i/,r) = e(^^+^^)>^(y)(l + o(l)), where ipp e Span{ijp : = /}, y^/j ^ 0. 
In the original variables (15. ip . the global patterns (I5.5p take the form: 

(5.6) Ufi{x,t) =t-^^ ip,3{jT7h)i^ + 0(1)) as t^+oo. 

5.2. Centre manifold patterns. Unlike the simpler case m = 1 in [20^, § 5.2], for the 
present m > 2, such patterns do exist. Performing a model ID analysis of the equation 
(15. 2p . as in [201 § 5], we conclude that such patterns may occur if 

(5.7) A/3 + ci = =^ / = = 2mci > 0, or p = pi. 
Studying the centre manifold behaviour of the simplest ID type 

(5.8) v{t) = ai(T)ipi + as r — )• +oo, 

we obtain from ( 15. 2 p the following equation for the expansion coefficient: 

(5.9) di = Ki\ai\P{l + o{l)), where m = , ^1) 0) , 

which admits global bounded orbits. For instance, noting that kq > 0, one obtains the 
behaviour 

(5.10) ao(r) = -[/«o(p- l)]"^^r"F^(l + o(l)) as r +oo {p = Po). 

Similarly, the same estimate is derived for any I > provided that ni ^ 0. Finally, 
this means that, at such critical values p = pi, we expect the following logarithmically 
perturbed patterns: 

I I JV+i N+l 

(5.11) uiix,t)^ -sign Ki[^] '-{tint) ^'"V^^(^) as t ^ +oo. 
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For the M- dimensional eigenspace for / > 1, we obtain the decomposition 

(5.12) v{y,T)= ap{T)'^p{y) + w^{y,T), 

that leads to a system of ODEs for the expansion coefficients {a^(r) }|/3|=;: 

(5.13) «7 = (I E|/3|=««/3(^)V^/3r,V^7> + ••• , |7|=^- 

Assuming the natural "homogenuity" of this centre subspace behaviour: 

(5.14) apir) = apT-^{l + o{l)) as r ^ +oo, \/3\=l, 

where {a^} are constants, f l5.13p reduces to an algebraic system (cf. the bifurcation one 
([229D) of the usual form: 

(5.15) ^, = -Mii\T.\^\=i^^M'^r,), \l\ = l- 

General solvability properties of f l5.15p . except some obvious elementary solutions, and 
sharp multiplicity results are unknown. Of course, as above, (15.151) is not variational. 

Again, as for m = 1, we arrive at two countable families of global patterns: the nonlinear 
(12. ip and the linearized (15.60 ones. Since the rescaled equation (15.20 is not a gradient 
system in any weighted space, their evolution completeness remains entirely open, though 
may be expected. 

References 

[1] R. Bari and B. Rynne, Solution curves and exact multiplicity results for 2mth order boundary value 

problems, J. Math. Anal. AppL, 292 (2004), 17-22. 
[2] M. Berger, Nonlinearity and Functional Analysis, Acad. Press, New York, 1977. 
[3] C.J. Budd, V.A. Galaktionov, and J.F. Williams, Self-similar blow-up in higher-order semilinear 

parabolic equations, SIAM J. Appl. Math., 64 (2004), 1775-1809. 
[4] M.S. Birman and M.Z. Solomjak, Spectral Theory of Self-Adjoint Operators in Hilbert Space, D. 

Reidel, Dordrecht/Tokyo, 1987. 
[5] R. Czaja and C. Rocha, Transversality in scalar reaction- diffusion equations on a circle, J. Differ. 

Equal., 245 (2008), 692-721. 
[6] E.A. Coddington and N. Levinson, Theory of Ordinary Differential Equations, McGraw-Hill Book 

Company, Inc., New York/London, 1955. 
[7] K. Deimling, Nonlinear Functional Analysis, Springer- Verlag, Berlin/Tokyo, 1985. 
[8] Yu.V. Egorov, V.A. Galaktionov, V.A. Kondratiev, and S.I. Pohozaev, On the necessary conditions 

of existence to a quasilinear inequality in the half-space, Comptes Rendus Acad. Sci. Paris, Serie I, 

330 (2000), 93-98. 

[9] Yu.V. Egorov, V.A. Galaktionov, V.A. Kondratiev, and S.I. Pohozaev, Global solutions of higher- 
order semilinear parabolic equations in the supercritical range, Adv. Differ. Equat., 9 (2004), 1009- 
1038. 

[10] S.D. Eidelman, Parabolic Systems, North-Holland Publ. Comp., Amsterdam/London, 1969. 
[11] U. Ellias, Eigenvalue problems for the equation Ly + Xp{x)y ~ 0, J. Differ. Equal., 29 (1978), 28-57. 

25 



[12] M. Escobedo and O. Kavian, Variational problems related to self-similar solutions of the heat equa- 
tion, Nonl. Anal., TMA, 11 (1987), 1103-1133. 

[13] B. Fiedler and C. Rocha, Connectivity and design of planar global attractors of Sturm type. II: 
Connection graphs, J. Differ. Equat., 244 (2008), 1255-1286. 

[14] M.V. Fedoryuk, Singularities of the kernels of Fourier integral operators and the asymptotic behaviour 
of the solution of the mixed problem, Russian Math. Surveys, 32 (1977), 67-120. 

[15] V.A. Galaktionov, Critical global asymptotics in higher-order semilinear parabolic equations. Int. J. 
Math. Math. Sci., 60 (2003), 3809-3825. 

[16] V.A. Galaktionov, Geometric Sturmian Theory of Nonlinear Parabolic Equations and Applications, 
Chapman &:HaU/CRC, Boca Raton, Florida, 2004. 

[17] V.A. Galaktionov and P.J. Harwin, Non-uniqueness and global similarity solutions for a higher-order 
semilinear parabolic equation, Nonlinearity, 18 (2005), 717-746. 

[18] V.A. Galaktionov, E. Mitidieri, and S.I. Pohozaev, Variational approach to complicated similar- 
ity solutions of higher-order nonlinear evolution equations of parabolic, hyperbolic, and nonlinear 
dispersion types. In: Sobolev Spaces in Mathematics. II, Appl. Anal, and Part. Differ. Equat., Se- 
ries: Int. Math. Scr., Vol. 9, V. Maz'ya Ed., Springer, New York, 2009 (an extended version in 
[arXrv:0902.1425 ). 

[19] V.A. Galaktionov, E. Mitidieri, and S.I. Pohozaev, Variational approach to complicated simi- 
larity solutions of higher-order nonlinear PDEs. II, Nonl. Anal.: RWA, 12 (2011), 2435-2466 
tarXiv: 1103.2643). 

[20] V.A. Galaktionov, E. Mitidieri, and S.I. Pohozaev, Classification of global and blow-up sign-changing 
solutions of a semilinear heat equation in the subcritical Fujita range I. Second-order diffusion, sub- 
mitted. 

[21] V.A. Galaktionov and S.I. Pohozaev, Existence and blow-up for higher-order semilinear parabolic 
equations: majorizing order-preserving operators, Indiana Univ. Math. J., 51 (2002), 1321-1338. 

[22] V.A. Galaktionov and J.L. Vazquez, A Stability Technique for Evolution Partial Differential Equa- 
tions. A Dynamical Systems Approach, Birkhauser, Boston/Berlin, 2004. 

[23] V.A. Galaktionov and J.F. Williams, On very singular similarity solutions of a higher-order semi- 
linear parabolic equation, Nonlinearity, 17 (2004), 1075-1099. 

[24] R. Joly and G. Raugel, Ceneric hyperbolicity of equilibria and periodic orbits of the parabolic equation 
on a circle. Trans. Amer. Math. Soc, 362 (2010), 5189-5211. 

[25] T. Kato, Perturbation Theory for Linear Operators, Springer- Verlag, Berlin/ New York, 1976. 

[26] M.A. Krasnosel'skii, Topological Methods in the Theory of Nonlinear Integral Equations, Pergamon 
Press, Oxford/Paris, 1964. 

[27] M.A. Krasnosel'skii and P.P. Zabreiko, Geometrical Methods of Nonlinear Analysis, Springer- Verlag, 
Berlin/Tokyo, 1984. 

[28] A. Lunardi, Analytic Semigroups and Optimal Regularity in Parabolic Problems, Birkhauser, 
Basel/Berlin, 1995. 

[29] B. Rynn, Global bifurcation for 2mth-order boundary value problems and infinitely many solutions 
of superlinear problems, J. Differ. Equat., 188 (2003), 461-472. 



26 



Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, UK 
E-mail address: vag@maths.bath.ac.uk 

Dipartimento di Matematica e Informatica, Universita di Trieste, Via Valerio 12, 34127 
Trieste, ITALY 

E-mail address: mitidier@units.it 

Steklov Mathematical Institute, Gubkina St. 8, 119991 Moscow, RUSSIA 
E-mail address: pokhozhaev@mi.ras.ru 



27 



